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Coherent  tunneling  between  two  InAs  quantum  dots  forms  delocalized  molecular  states.  Using 
magnetophotoluminescence  spectroscopy  we  show  that  when  holes  tunnel  through  a  thin  barrier,  the 
lowest  energy  molecular  state  has  bonding  orbital  character.  However,  as  the  thickness  of  the  barrier 
increases,  the  molecular  ground  state  changes  character  from  a  bonding  orbital  to  an  antibonding  orbital, 
confirming  recent  theoretical  predictions.  We  explain  how  the  spin-orbit  interaction  causes  this  counter¬ 
intuitive  reversal  by  using  a  four-band  k  ■  p  model  and  atomistic  calculations  that  account  for  strain. 
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Quantum  dots  have  confined  energy  levels  analogous  to 
ordinary  atoms.  Two  quantum  dots  in  close  proximity  can 
be  viewed  as  an  artificial  diatomic  molecule  when  coherent 
tunnel  coupling  leads  to  the  formation  of  delocalized 
states.  The  properties  of  such  quantum-dot  molecules 
(QDMs)  have  been  the  focus  of  much  research  because 
of  potential  applications  in  novel  optoelectronic  devices  or 
quantum  information  processing.  In  analogy  with  natural 
diatomic  molecules,  one  expects  the  lowest  energy  delo¬ 
calized  molecular  state  to  have  bonding  orbital  character. 
However,  recent  theoretical  studies  have  predicted  that  the 
molecular  ground  state  for  a  hole  in  an  InAs  QDM  can  have 
antibonding  character  [1-3].  If  verified  by  experiment,  an 
antibonding  molecular  ground  state  would  provide  a  strik¬ 
ing  example  of  a  novel  property  of  artificial  atoms  that 
cannot  simply  be  explained  as  a  rescaled  version  of  the 
physics  of  real  atoms. 

In  this  Letter  we  present  the  first  experimental  observa¬ 
tion  of  an  antibonding  molecular  ground  state.  We  find  that 
the  molecular  ground  state  changes  character  from  a  bond¬ 
ing  orbital  to  an  antibonding  orbital  as  the  thickness  of  the 
barrier  separating  the  two  coupled  quantum  dots  is  in¬ 
creased.  Using  a  four-band  k  ■  p  model  validated  by  atom¬ 
istic  calculations,  we  explain  how  this  counterintuitive 
result  arises  from  the  spin-orbit  (SO)  interaction. 

We  use  magneto-optical  spectroscopy  to  study  QDMs 
composed  of  two  vertically  stacked  InAs/GaAs  quantum 
dots  separated  by  a  GaAs  barrier.  The  two  dots  have  differ¬ 
ent  size,  composition,  and  strain,  and  therefore  different 
confined  energy  levels.  As  a  result,  the  electron  and  hole 
tend  to  localize  in  individual  dots,  as  depicted  in  the  left- 
hand  insets  of  Fig.  1.  Delocalized  molecular  states  are 
formed  by  coherent  tunneling  [4]  when  an  electric  field 
tunes  the  relative  energies  of  confined  states  in  the  two  dots 
through  resonance  [5-7].  Either  electron  or  hole  tunneling 
can  be  induced  [8],  but  in  this  work  we  focus  only  on  hole 
tunneling.  Because  of  the  large  inhomogeneous  distribu¬ 
tion  of  parameters  in  ensembles  of  QDMs,  all  spectroscopy 
is  performed  on  single  QDMs. 


PACS  numbers:  78.67.Hc,  78.20.Ls,  78.47.-p,  78.55.Cr 

Figure  1  shows  the  anticrossing  of  the  neutral  exciton 
(Y°)  that  results  from  coherent  tunneling  of  a  single  hole 
through  a  thin  (2  nm)  barrier  while  the  electron  remains 
localized  in  the  bottom  dot.  The  tunneling  of  holes  creates 
molecular  states  that  are  the  symmetric  and  antisymmetric 
combinations  of  the  two  basis  states  where  the  hole  is  in 
one  dot  or  the  other  [9].  In  analogy  to  real  molecules,  we 
call  the  symmetric  (nodeless)  molecular  state,  which  has  an 
enhanced  wave  function  amplitude  in  the  barrier,  a  bonding 
state.  The  antisymmetric  (noded)  state  has  a  suppressed 
amplitude  in  the  barrier  and  is  called  the  antibonding  state. 
These  molecular  orbitals  are  depicted  schematically  by  the 
right-hand  insets  in  Fig.  1.  Intuitively  one  expects  the 
molecular  ground  state  to  have  bonding  orbital  character 
and  the  first  excited  molecular  state  to  have  antibonding 
orbital  character. 

The  formation  of  molecular  orbitals  at  an  anticrossing  is 
described  by  a  simple  Hamiltonian  using  an  atomiclike 


FIG.  1  (color  online).  Photoluminescence  (PL)  measurement 
of  the  electric  field-induced  anticrossing  of  X°  at  zero  magnetic 
field  for  a  sample  with  2  nm  barrier.  A  indicates  the  anticrossing 
energy  gap.  Insets:  If  the  hole  energy  levels  are  out  of  resonance 
(left)  the  hole  is  localized  in  one  of  the  individual  dots.  When  the 
hole  levels  are  tuned  into  resonance  by  the  applied  electric  field, 
coherent  tunneling  leads  to  the  formation  of  bonding  (bottom 
right)  and  antibonding  (upper  right)  molecular  wave  functions. 
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basis  with  a  hole  either  in  one  dot  or  the  other: 

M-.  *-"/'+/„}  (|) 

Here  E0  is  the  energy  of  the  localized  hole  states  at 
resonance,  t  is  the  tunneling  rate,  and  /  =  edF  is  the 
Stark  energy  due  to  the  electric  field  F.  The  energies  of 
the  bonding  and  antibonding  molecular  orbitals  are  given 
by  the  eigenvalues  of  Eq.  (1).  When  the  electric  field  is 
tuned  to  resonance  (/  =  /0),  the  energies  are  Eh  =  E0  —  t 
and  Eah  =  E0  +  t.  The  magnitude  of  t  is  determined  by  the 
splitting  between  the  two  molecular  states:  2t  =  Eab  —  Eb. 
The  sign  of  t  is  determined  by  which  orbital  state  is  at 
higher  energy.  However,  the  sign  of  t  cannot  be  measured 
from  anticrossing  energy  gaps  (A  =  2\t\)  in  photolumines¬ 
cence  (PL)  spectra  like  that  of  Fig.  1 . 

Using  magneto-PL  we  can  directly  measure  the  orbital 
character  of  the  molecular  states  and  determine  the  sign  of 
t  [10].  When  a  magnetic  field  is  applied  to  InAs/GaAs 
QDMs  with  a  2  nm  GaAs  barrier  there  is  a  large  resonant 
change  in  the  Zeeman  energy  splitting  as  a  function  of 
electric  field  [10].  This  effect  is  shown  in  Fig.  2(a),  where 
we  plot  the  PL  lines  of  the  same  QDM  shown  in  Fig.  1, 
but  now  in  a  longitudinal  magnetic  field  of  B  =  6  T.  The 
resonant  change  in  Zeeman  splitting  is  plotted  in  Fig.  2(c). 
The  black  [gray  (red)]  shading  indicates  the  resonant 
change  in  Zeeman  splitting  for  the  molecular  ground  (ex¬ 
cited)  state.  As  discussed  below,  the  resonant  changes  in 


b 


Electric  Field  (kV/cm) 


FIG.  2  (color  online).  Zeeman  energy  splitting  as  a  function  of 
applied  electric  field  for  B  =  6  T.  (a),(b)  Energies  of  the  X°  PL 
lines  for  QDMs  with  2  and  4  nm  barriers.  Solid  lines  and  dashed 
lines  calculated  using  Eq.  (2)  indicate  the  two  separate  spin 
configurations,  (c)-(f)  QDMs  have  a  barrier  thickness  of  (c)  2, 
(d)  3,  (e)  4,  and  (f)  6  nm.  Solid  curves  are  calculated  with  Eq.  (3) 
using  (h0,  ti)  =  (-0.368,  0.229),  (-0.446,  0.082),  (-0.404, 
0.076)  meV  for  2,  3,  and  4  nm,  respectively. 


Zeeman  splitting  arise  from  the  contribution  of  the  GaAs 
barrier  to  the  net  g  factor  for  the  delocalized  hole  [10]. 

In  Figs.  2(d)-2(f)  we  show  the  measured  change  in 
Zeeman  splitting  from  samples  with  increasing  barrier 
thickness  d.  When  d  =  3  and  4  nm  [Figs.  2(d)  and  2(e)], 
there  is  a  dramatic  reversal  in  the  nature  of  the  resonance: 
the  molecular  ground  state  (black  shading)  now  shows  a 
resonant  increase  in  Zeeman  splitting  and  the  first  molecu¬ 
lar  excited  state  [gray  (red)  shading]  now  shows  the  reso¬ 
nant  decrease  in  Zeeman  splitting.  The  amplitude  of  the 
resonant  change  in  Zeeman  splitting  decreases  as  the 
thickness  of  the  barrier  increases,  and  is  below  our  noise 
level  for  d  =  6  nm  [Fig.  2(f)].  The  decreasing  amplitude 
results  from  the  reduction  of  the  amplitude  of  the  wave 
function  in  the  barrier  with  increasing  barrier  thickness. 
We  will  first  show  from  an  analysis  of  the  data  that  the 
inversion  of  the  Zeeman  resonance  can  be  understood  as  a 
change  in  the  sign  of  t.  Then  we  will  discuss  how  this  can 
be  understood  from  theoretical  considerations. 

To  quantitatively  analyze  the  resonant  changes  in 
Zeeman  splitting  (Fig.  2)  and  the  associated  tunneling 
rates,  we  add  the  Zeeman  interaction  to  Eq.  (1): 
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We  obtain  two  Hamiltonians  for  the  two  spin  configura¬ 
tions  ]]]  and  |]L ,  where  (|,  |)  and  (]],  ft)  are  the  electron  spin 
and  hole  spinor  projections,  respectively.  The  diagonal 
part,  h0  =  ( ge  +  gh)/jLBB/2,  is  just  the  normal  Zeeman 
energy  for  the  isolated  dots.  ge  and  gh  are  the  g  factors 
for  electrons  and  holes.  The  off-diagonal  part,  /?/,  creates  a 
spin-dependent  contribution  to  the  tunneling  rate  (—  t  => 
—  t  ±  h'),  which  is  responsible  for  the  resonant  change  in 
Zeeman  splitting  [11]. 

To  obtain  an  expression  for  the  Zeeman  splitting  of  the 
ground  (Ac)  and  excited  (A£)  states,  we  take  the  difference 
between  the  eigenvalues  for  the  two  spin  states  in  Eq.  (2): 
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Using  Eq.  (3)  we  obtain  the  fits  to  the  Zeeman  resonance 
data  shown  in  Figs.  2(c)-2(f).  The  off-resonant  value  of  the 
splitting  is  determined  by  h0  and  the  magnitude  of  the 
resonance  is  given  by  h' .  The  fit  values  are  given  in  the 
figure  caption.  The  width  of  the  resonance  is  determined  by 
the  magnitude  of  the  tunneling  rate  t,  which  is  indepen¬ 
dently  determined  by  the  anticrossing  energy  (for  example, 
A  in  Fig.  1). 

This  analysis  leads  to  a  simple  intuitive  description  of 
the  origin  of  the  resonant  change  in  Zeeman  splitting.  As 
the  electric  field  tunes  the  states  of  the  two  dots  into 
resonance,  the  formation  of  molecular  orbitals  with  bond¬ 
ing  or  antibonding  character  changes  the  amplitude  of  the 
molecular  hole  wave  function  in  the  barrier,  which  has  a 
different  hole  g  factor  [10],  The  second  term  on  the  right- 
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hand  side  of  Eq.  (3)  captures  the  resonant  change  in  the 
contribution  from  the  barrier  (h1)  and  is  responsible  for  the 
resonant  change  in  the  Zeeman  splitting.  The  relative  signs 
of  the  first  and  second  terms  in  Eq.  (3)  determine  whether 
there  is  a  resonant  enhancement  or  suppression  of  the 
splitting.  We  find  that  h0  is  negative  [12].  Thus  if  h't  is 
positive,  the  Zeeman  splitting  of  the  ground  state  is  sup¬ 
pressed  at  resonance  and  the  splitting  of  the  excited  state  is 
enhanced.  This  is  what  we  measure  for  the  d  —  2  nm 
sample.  We  will  show  below  from  theory  that  t  >  0  for 
thin  barriers,  corresponding  to  the  normal  case  where  the 
ground  state  is  a  bonding  orbital.  Therefore  h'  is  also 
positive. 

We  make  the  assumption  that  the  sign  of  h'  does  not 
change  as  the  bander  thickness  is  increased.  This  assump¬ 
tion  is  reasonable:  as  the  banier  becomes  thicker,  the  sign 
of  h'  should  tend  toward  the  g  factor  for  holes  in  bulk 
GaAs,  which  is  known  to  be  positive  [13,14],  Because  h' 
remains  positive,  it  is  a  change  of  sign  of  t  in  Eq.  (3)  that 
leads  to  the  inversion  of  the  resonant  Zeeman  energy 
splitting  as  the  barrier  thickness  is  increased.  This  reversal 
in  the  sign  of  t  means  that  there  is  a  reversal  of  the  energy 
of  the  bonding  and  antibonding  states,  i.e.,  the  antibonding 
state  becomes  the  energetic  ground  state.  We  now  consider 
the  theoretical  origin  of  such  a  reversal  and  show  that  it  is  a 
result  of  the  SO  interaction. 

Holes  experience  a  strong  SO  interaction  because  they 
are  derived  from  /j-typc  atomic  orbitals  of  the  semicon¬ 
ductor  lattice  [15].  In  QDMs,  the  SO  interaction  couples 
the  hole’s  atomic  orbital  and  spin  degrees  of  freedom  to 
give  a  total  (Bloch)  angular  momentum  J  =  3/2.  Jz  = 
±3/2  projections  correspond  to  heavy  holes  and  Jz  = 
±1/2  to  light  holes.  Because  the  light-hole  states  are 
shifted  up  in  energy  by  confinement  and  strain,  it  is  often 
useful  to  represent  the  two  low-energy  (heavy-hole)  states 
as  pseudo-spin- 1/2  particles  (]],!]).  When  the  heavy-light- 
hole  mixing  is  included,  low-energy  hole  states  can  be 
described  within  the  Luttinger-Kohn  k  ■  p  Hamiltonian 
formalism  as  four  component  Luttinger  spinors  [3].  Each 
spinor  is  an  admixture  of  all  four  projections  of  Jz,  with  the 
heavy  hole  typically  the  dominant  component.  However,  as 
we  show  here,  minor  components  of  no  more  than  5%  are 
sufficient  to  substantially  alter  the  character  of  the  molecu¬ 
lar  orbitals. 

The  influence  of  the  minor  components  is  apparent  from 
Fig.  3.  In  Fig.  3(a),  we  plot  the  energies  of  the  molecular 
ground  and  first  excited  states  for  asymmetric  QDMs  sub¬ 
ject  to  a  resonant  electric  field  calculated  using  a  simple 
one-band  effective  mass  model,  which  neglects  SO  inter¬ 
actions.  As  expected,  the  energy  separation  of  the  bonding 
(solid  black  line)  and  antibonding  (dashed  red  fine)  states 
decreases  as  a  function  of  increasing  barrier  thickness  and 
the  bonding  orbital  remains  the  molecular  ground  state.  In 
Fig.  3(b)  we  show  the  energies  of  the  bonding  (solid  blue 
line)  and  antibonding  (dashed  red  line)  states  calculated 


— - - - -  ■  ,  ■  -1-1— T— - , - , - , - , - 

01234567  01234567 


Barrier  Thickness  d  (nm)  Barrier  Thickness  d  (nm) 

FIG.  3  (color  online).  (a),(b)  Energy  of  the  ground  and  first 
excited  molecular  states  as  a  function  of  d  calculated  using 
(a)  the  single-band  effective  mass  and  (b)  k  •  p  theory.  Scale 
bars  are  5  meV.  Insets  show  the  orbital  character  of  the  dominant 
hole  spinor  component,  (c)  Tunneling  rates  of  a  single  hole 
versus  d  calculated  as  described  in  the  text.  Inset:  Schematic 
depiction  of  a  QDM.  (d)  Schematic  depiction  of  the  SO  induced 
mixing  between  bonding  ( Eh )  and  antibonding  (Ea)  molecular 
states  of  the  heavy  (EH)  and  light  ( EL )  holes. 


using  a  four-band  k  ■  p  model  that  includes  the  SO  inter¬ 
action.  At  d  ~  1.75  nm  the  energies  of  the  bonding  and 
antibonding  states  cross  and  the  antibonding  state  becomes 
the  molecular  ground  state.  We  estimate  the  antibonding 
character  of  the  ground  state  spinor  for  large  barrier  thick¬ 
nesses  to  be  as  large  as  95%,  many  times  larger  than  in 
known  atomic  systems  [16]. 

In  Fig.  3(c)  we  plot  the  values  for  t0  and  t,  i.e.,  half  the 
difference  between  the  state  energies  given  in  Figs.  3(a) 
and  3(b),  respectively.  In  the  absence  of  SO  interaction,  the 
tunneling  rate  is  determined  simply  by  the  overlap  between 
the  hole  orbitals  of  the  individual  dots  (t0),  which  decreases 
exponentially  with  increasing  barrier  thickness  at  a  rate 
dependent  on  the  heavy-hole  mass.  When  the  SO  interac¬ 
tion  is  included,  there  is  a  correction  to  the  tunneling  rate, 
t  =  t0  —  fSo-  This  tso  term  arises  from  the  small  contribu¬ 
tion  of  the  light-hole  component  of  the  spinor.  The  light- 
hole  component  has  approximate  parity  along  z  opposite  to 
that  of  the  heavy-hole  component  [3,16].  The  light-hole 
component  therefore  adds  a  small  antibonding  (bonding) 
component  to  the  bonding  (antibonding)  state  determined 
by  the  dominant  heavy-hole  component,  as  shown  sche¬ 
matically  in  Fig.  3(d).  The  addition  of  this  antibonding 
component  increases  the  energy  of  the  bonding  state  and 
vice  versa  for  the  antibonding  state.  As  the  barrier  thick¬ 
ness  increases,  tso  does  not  decrease  as  fast  as  t0,  in  part 
because  of  its  light-hole  origin.  For  thin  barriers,  the  tso 
correction  is  small  compared  to  the  large  r(l,  and  t  remains 
positive.  However,  for  thicker  barriers  t0  decreases  and 
becomes  comparable  to  tso.  When  t0  <  tso  the  tunneling 
rate  is  negative  and  the  antibonding  orbital  is  the  molecular 
ground  state. 

We  have  verified  that  this  simple  four-band  k  ■  p  ap¬ 
proach  captures  the  essential  physics  of  the  system  by 
comparison  with  an  atomistic  calculation  (using  >106 
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FIG.  4  (color  online),  (a)  Molecular  ground  states  for  each 
charge  configuration  are  determined  by  sequentially  filling  mo¬ 
lecular  orbitals  with  holes,  (b)  Calculated  energy  levels  for  the 
X2+  transition  (from  three  holes  plus  one  electron  to  two  holes) 
and  the  X°  transition  (from  one  hole  plus  one  electron  to  zero 
holes),  (c)  Electric  field  dependence  of  the  Zeeman  splitting  for 
the  X°  and  X2+  molecular  ground  states  in  a  QDM  with  2  nm 
barrier  at  B  =  6  T.  The  resonances  peak  at  two  different  values 
of  the  electric  field  (fx o  and  fx 2+)  because  of  different  Coulomb 
interactions. 


tal  character.  The  three-hole  state  is  observed  as  the  initial 
state  of  the  doubly  charged  exciton  (three  holes  and  one 
electron)  [16,17], 

Here  we  have  presented  the  first  experimental  ob¬ 
servation  of  an  antibonding  molecular  ground  state. 
Antibonding  molecular  ground  states  are  never  observed 
in  natural  molecules,  so  this  result  provides  a  striking 
example  of  the  new  properties  that  can  be  engineered  using 
semiconductor  nanostructures.  Our  result  specifically  dem¬ 
onstrates  that  the  SO  interaction  can  be  used  to  design  hole 
molecular  ground  states  with  arbitrary  orbital  character. 
Possible  applications  in  spintronics  and  quantum  informa¬ 
tion  processing  include  designing  structures  and  protocols 
to  coherently  control  single  confined  spins  with  electric 
fields  or  g  tensor  modulation  [18]. 
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atoms)  of  the  hole  levels  of  a  QDM  described  by  the 
sp3d5s*  tight-binding  model  [16].  This  approach  accounts 
for  strain  and  changes  to  the  underlying  crystal  lattice  on 
the  atomistic  level.  The  results  of  this  calculation  [solid 
blue  points  in  Fig.  3(c)]  qualitatively  match  the  k  ■  p 
results. 

Our  theoretical  model  predicts  an  antibonding  molecu¬ 
lar  ground  state  for  a  barrier  thickness  d  >  1.75  nm. 
Experimentally,  we  find  that  all  examples  (7)  in  the  sample 
with  d  =  2  nm  have  a  molecular  ground  state  with  bond¬ 
ing  orbital  character;  in  the  sample  with  d  =  4  nm,  all 
examples  (3)  have  an  antibonding  ground  state.  The  inter¬ 
mediate  case  {d  =  3  nm)  has  examples  of  both  types  of 
behavior,  indicating  that  the  reversal  of  orbital  character 
occurs  near  <7=3  nm.  The  small  discrepancy  with  theory 
most  likely  arises  from  details  of  dot  structure,  and  the 
coexistence  of  both  behaviors  at  d  =  3  nm  most  likely 
arises  from  fluctuations. 

We  now  consider  what  happens  when  we  add  additional 
holes  to  the  molecular  orbitals.  In  Fig.  4(a)  we  schemati¬ 
cally  depict  the  filling  of  the  molecular  orbitals  when  the 
bonding  state  is  the  lowest  energy  single  particle  state  (i.e., 
d  =  2  nm).  The  lowest  energy  two-hole  state  is  a  spin 
singlet  with  bonding  character  that  has  no  magnetic  field 
spin  splitting.  The  lowest  energy  three-hole  state  must  have 
an  unpaired  hole  in  the  antibonding  orbital.  Because  of  the 
unpaired  hole,  this  state  should  have  a  magnetic  field 
splitting  like  that  of  the  one-hole  state,  but  the  Zeeman 
splitting  should  increase  on  resonance  because  the  un¬ 
paired  hole  is  in  the  antibonding  orbital.  This  is  exactly 
what  we  observe,  as  shown  in  Fig.  4(c).  The  one-hole  and 
three-hole  states  of  the  same  QDM  (with  d  =  2  nm)  have 
opposite  behavior,  indicating  that  they  have  different  orbi- 
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